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EXECUTIVE  SUMMARY 


The  1989  JASON  Summer  Study  on  Drag  Reduction  focused  on  the 
physics  which  underlies  methods  utilizing  polymer  additives.  The  study  in¬ 
cluded  a  review  of  drag  reduction  phenomenology,  the  development  of  con¬ 
tinuum  models  of  the  dynamics  of  dilute  polymer  solutions,  the  introduction 
of  a  simple  theory  of  polymer-hydrodynamic  turbulence  and  a  discussion 
of  its  implications  for  grid  flows,  and  consideration  of  optimizing  polymer 
architecture. 

Current  experimental  evidence  from  studies  of  wall  streak  spacing  and 
wall  stress  attenuation  suggests  that  drag  reduction  occurs  by  the  quench¬ 
ing  of  small  scale,  high  frequency  eddys  in  a  buffer  layer  adjacent  to  the 
viscous  sub-layer.  Theoretical  interpretations  purport  to  explain  drag  reduc¬ 
tion  in  terms  of  a  buffer  layer  viscosity  enhanced  by  flow-induced  polymer 
uncoiling(G.  Ryskin),  or  in  terms  of  a  truncated  viscoelastic  cascade  (P.G. 
DeGennes  and  M.  Tabor).  The  uncoiling  model  is  open  to  question,  because 
it  is  based  on  the  application  of  intuition  from  static  extensional  flows  to 
fully  developed  turbulence.  The  viscoelastic  model  is  incomplete  because  it 
fails  to  consistently  account  for  the  ultimate  disposition  of  turbulent  eddy 
energy. 

As  a  foundation  for  further  theoretical  and  computational  study,  a  set  of 
macroscopic  continuum  equations  which  self-consistently  describe  the  fluid 
and  polymer  dynamics  have  been  derived.  In  their  simplest  form,  appro¬ 
priate  to  the  linear  “dumbell”  polymer  model,  these  equations  describe  the 
dynamics  of  the  polymer  extension  field  and  the  fluid  velocity  field,  including 
back-reaction  effects  due  to  the  presence  of  the  polymers.  The  continuum 
model  is  quite  similar  to  magnetohydrodynamics,  with  the  polymer  exten¬ 
sion  field  analogous  to  the  magnetic  field.  In  addition  to  dissipation  due 
to  viscosity  and  polymer  diffusion,  a  scale-independent  damping  term,  rep¬ 
resenting  the  effects  of  polymer  relaxation  modes,  appears  in  the  extension 
field  equation.  We  have  found  the  analogy  between  polymer-hydrodynamic 
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turbulence  and  magnetohydrodynamic  turbulence  quite  useful.  The  exten¬ 
sion  field  is  excited  by  fluid  eddy  stretching,  until  elastic  energy  approaches 
fluid  kinetic  energy.  This  stretching  process  is  limited  by  polymer  relaxation 
modes  at  low  frequencies. 

When  excited,  the  extensional  field  inhibits  the  eddy  interaction,  resulting 
in  a  change  of  the  inertial  range  spectrum  from  k~^  to  k~^.  Direct  calcula¬ 
tions  then  indicate  the  decay  rate  of  grid  polymer-hydrodynamic  turbulence 
exceeds  that  of  ordinary  fluid  grid  turbulence.  This  is  due  (primarily)  to  dis¬ 
sipation  due  to  polymer  relaxation  modes.  Extrapolation  to  the  case  of  wall 
flows  suggests  a  picture  of  a  drag-reduced  flow  as  a  Newtonian  slug  with  a 
visco-elastic  buffer  layer.  Enhanced  dissipation  is  due  to  polymer  relaxation 
modes. 

Consideration  was  given  to  optimizing  polymer  architecture  for  the  pur¬ 
pose  of  drag  reduction.  While  linear  chain  polymers  best  meet  the  dual  cri¬ 
teria  of  low  relaxation  frequency  and  high  permeated  fluid  volume  fraction 
(by  weight),  branched  polymers  can  better  resist  degradation  by  distributing 
(mid-coil)  flow  induced  stresses. 
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1  OVERVIEW 


Polymer  drag  reduction  was  discovered  about  forty  years  ago  by  Toms[l], 
who  observed  drag  reductions  of  30  —  40  percent  upon  adding  only  lOppm 
by  weight  of  polymer(methylmethacrylate)  to  turbulent  monochlorobeneze 
flowing  down  a  pipe.  Similar  drag  reductions  have  been  observed  in  poly¬ 
isobutylene  polymer  in  benzene  and  cyclohexane.  For  many  applications,  the 
most  important  solvent  is  undoubtedly  water,  for  which  polyethylene  oxide 
(PEO)  produces  very  similar  effects.  Minute  concentrations  of  PEO  have 
been  used  to  increase  the  capacities  of  irrigation  networks,  and  municipal 
sewer  systems  and  to  extend  the  range  of  fire  hoses.  The  drag  reducing 
effects  of  these  polymers  are  most  pronounced  when  the  polymers  are  flex¬ 
ible  with  high  molecular  weights,  corresponding  to  iV  =  10^  or  N  =  10® 
individual  monomer  units.  Despite  the  great  technological  importance  of 
these  materials,  there  is  as  yet  no  fundamental  understanding  of  how  this 
remarkable  drag  reduction  is  achieved.  In  this  section  we  briefly  discuss  the 
experiments  and  some  important  elementary  facts  about  polymer  sizes  and 
relaxation  rates.  Some  novel  chemical  architectures  for  polymer  additives 
are  discussed.  Section  2  is  a  critique  of  existing  ideas  about  polymer  drag 
reduction.  In  Section  3.1  we  discuss  macroscopic  continuum  models  of  dilute 
solutions  of  polymers  modeled  as  stretchable  dumbells.  In  Section  3.2,  we 
discuss  the  dynamics  of  polymer-hydrodynamic  (PHD)  turbulence.  The  con¬ 
tinuum  model  presented  in  Section  3  can  be  analyzed  using  methods  fzuniliar 
from  studies  of  turbulent  magneto-hydrodynamics.  The  analysis  suggests  in¬ 
teresting  and  non-trivial  modifications  of  the  Kolmogorov  turbulence  theory. 
A  prediction  of  the  energy  dissipation  rate  for  grid  PHD  turbulence  is  given, 
and  its  implications  for  pipe  flow  discussed.  Section  4  contains  our  program 
recommendations.  The  basic  phenomenology  of  MHD  turbulence  is  reviewed 
in  the  Appendix. 
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1.1  Phenomenology  of  Drag  Reduction 


A  typical  experimental  result  [2]  is  shown  in  Figure  1-1.  For  a  fluid  with 
density  p  and  average  velocity  <v>  flowing  down  a  pipe  of  diameter  D  the 
figure  of  merit  is  the  mass  flux 

Q  =  <v>  (1-1) 
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obtained  for  a  given  pressure  head  Ap.  A  dimensionless  measure  of  the  ratio 
of  pressure  head  to  mass  flux  is  the  “friction  factor” 


DAp 

2p  <v>^ 


(1-2) 


which  is  /  =  16/i2e  for  laminar  flows  with  Reynolds  number 
Re  =  D  <v  >  jp.  Although  the  addition  of  small  amounts  of  PEO  polymer 
has  no  observable  effect  in  the  laminar  regime,  there  is  a  striking  change 
in  the  drag  above  the  transition  to  turbulence  which  occurs  for  Re  «  2500 
obtainable  with  as  little  as  5ppm  (by  weight)  of  additive.  The  effect  increases 
with  increasing  polymer  concentration. 


It  is  well  known  that  enhanced  drag  in  turbulent,  high  Reynolds  number 
flows  is  caused  by  radial  transport  of  flow  momentum  by  fluid  eddys.  Stan¬ 
dard  pipe  flows  are  well  described  by  the  Prandtl  mixing  length  theory[4], 
in  which  the  flow  is  sub-divided  into  a  viscous  sublayer  near  the  wall  and 
an  inertial  layer  extending  to  the  core  of  the  pipe.  In  the  viscous  sublayer 
(0  <  X  <  vfU^j  where  v  is  the  fluid  kinetic  viscosity  and  C/,  the  friction  veloc¬ 
ity),  the  flow  profile  is  linear  {U  =  U^xlv)  and  all  eddys  fall  in  the  dissipation 
range.  In  the  inertial  layer  {u/U,  <  x),  the  flow  profile  follows  the  logarith¬ 
mic  “Law  of  the  Wall”  U  =  {U»fK)(nx  {k  is  the  von  Karman  constant),  with 
a  spectrum  of  turbulent  eddys  of  size  f  <  x  to  be  found  at 

each  distance  x  from  the  wall.  At  the  onset  of  polymer  drag  reduction,  an 
elastic  buffer  layer  adjacent  to  the  viscous  sublayer  arises.  Within  the  buffer 
layer,  the  flow  profile  steepens  relative  to  the  Newtonian  logarithmic  profile, 
resulting  in  the  formation  of  a  region  of  “effective  slip”  (i.e.  reduced  momen¬ 
tum  transport).  It  is  interesting  to  note  that  in  the  buffer  layer,  fluctuating 
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'  Friction  factor  f 


Figure  1-1.  Plot  of  the  "friction  factor"  defined  in  the  text  as  a  function  of  Reynolds  number  for 
turbulent  water  flowing  in  a  pipe  with  various  concentrations  of  PEO.  Points 
labelled  "N"  are  with  no  polymer  additives,  while  points  labelled  "P”  have 
polymers  in  varying  concentrations;  from  Ref.  [2] 
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radial  velocities  are  smaller  than  their  Newtonian  counterparts,  but  fluctu¬ 
ating  axial  velocities  actually  increase!  However,  radial  and  axial  velocities 
tend  to  decouple  in  drag  reduced  flows,  so  that  <  VU  >  decreases  and  the 
net  turbulent  transport  of  momentum  is  reduced[3].  Note  that  the  increase 
in  axial  fluctuation  energy  is  consistent  with  the  formation  of  a  region  of  en¬ 
hanced  slip.  As  drag  reduction  increases,  the  buffer  layer  expands  tow^Lrd  the 
center  of  the  channel  while  the  inertial  layer  shrinks.  It  is  important  to  keep 
in  mind  that  the  viscous  sublayer  remains  intact  throughout  the  drag  reduc¬ 
tion  process.  Thus,  drag  reduced  flows  can  be  viewed  as  a  Newtonian  plug, 
connected  to  the  viscous  sublayer  by  an  elastic  buffer  layer  which  supports 
enhanced  slip. 

It  is  worthwhile  to  comment  on  the  observations  pertaining  to  the  onset  of 
drag  reduction.  Experimental  studies  indicate  that  a  critical  wall  shear  stress 

is  required  for  the  onset  of  drag  reduction,  and  that  <7^,  is  independent  of 
pipe  diameter,  solvent  density,  or  viscosity.  A  dependence  of  <7^,  on  polymer 
gyration  radius  (defined  in  Section  1.2)  Rq  of  the  form  cr^  ^  Rq^  is  observed, 
with  2  <  a  <  3.  The  gyration  radius  Ra  is  virtually  always  significantly 
smaller  than  the  dissipation  scales  of  turbulent  flows  of  Reynolds  numbers 
Re  <  10®. 

There  are  also  experiments[5]  which  monitor  the  friction  factor  as  a  func¬ 
tion  of  the  distance  downstream  from  a  point  in  the  center  of  the  pipe  where 
the  drag  reducing  agent  is  injected.  Some  drag  reduction  appears  even  be¬ 
fore  significant  amounts  of  additive  have  had  time  to  diffuse  to  the  wall.  One 
(controversial)  interpretation  of  these  experiments  is  that  polymers  do  not 
merely  lubricate  the  boundary  lay^  at  the  wall  but  also  act  to  modify  the 
approximately  homogeneous  isotropic  turbulence  in  the  center  of  the  pipe. 
We  are  unaware  of  experiments  on  changes  induced  by  polymers  added  to 
grid-generated  turbulence,  away  from  walls,  which  could  clarify  this  point. 
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1.2  Elementary  Facts  About  Polymers 


To  get  some  feeling  for  the  ways  in  which  polymer  additives  might  affect 
turbulent  flow,  we  need  to  estimate  their  characteristic  size  and  time  scale[6]. 
A  simple  random  walk  model  of  a  polymer  chain  with  N  monomer  units, 
which  ignores  the  interactions  between  different  monomers,  predicts  a  size  or 
“radius  of  gyration”  Rg  which  grows  like  A?, 

RgfiiaNK  (1-3) 

where  a  is  the  monomer  size.  A  more  accurate  formula,  which  takes  into 
account  swelling  due  to  self-avoiding  interactions  between  monomers  is, 

RgfuaNK  (1-4) 


The  chemical  formula  for  polyethylene  oxide  with  N  monomers  is 

. .  .\—CH2  CH2  ~  0—]ff . . . .  (1—5) 

Each  monomer  unit  is  about  a  «  5A.  The  density  of  polymeric  materiaJ  with 
iV  =  10^  inside  a  cube  with  edge  Rg  centered  on  this  very  tenuous  object  is 
thus  only  4  x  10~*g/cm^l.  This  remarkably  low  density  of  monomer  means 
that  no  water  molecule  is  ever  very  far  from  a  flexible  polymer  chain  even 
at  very  low  concentrations  of  chains.  The  volume,  for  example,  enclosed  by 
the  radii  of  gyration  of  1.4  x  10*^  PEO  molecules  with  N  =  10^  added  to 
one  cubic  centimeter  of  water  (corresponding  to  a  concentration  of  lOOppm 
by  weight)  is  about  25  percent  of  the  available  volume. 


Each  polymer  ch^un  is  a  delicate  little  mechanical  transponder,  with  its 
own  set  of  internal  relaxation  times,  which  reacts  to  the  stresses  associated 
with  fluid  motion.  Although  the  spectrum  of  relaxation  times  of  a  polymer 
in  a  solvent  is  complicated,  the  slowest  dynamical  mode  is  associated  with 
distortions  of  the  (approximately  spherical)  equilibrium  polymer  shape  and 
is  known  to  be[6] 


Tz 


ksT 


^1.8 


(1-6) 
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where  j/,  is  the  solvent  viscosity.  The  important  time  scale  r,  is  called  the 
“Zimm  time”  and  is  of  order  milliseconds  for  PEO  v/ith  N  —  10^  monomer 
units.  To  within  numerical  factors,  it  can  be  written  as[6] 

r,  «  'ilK  (1  -  7) 

i.e.,  as  the  ratio  of  the  Stokes  drag  coefficient  7  =  &‘Kr},Rg  for  a  slightly 
stretched  polymer  of  characteristic  size  Rg  and  a  Hookean  force  constant 
K  =  ksT /Rg.  Temperature  appears  in  the  effective  polymer  spring  constant 
because  the  restoring  force  for  a  stretched  polymer  is  entropic  in  origin,  due  to 
the  many  more  configurations  available  to  a  polymer  in  its  coiled,  as  opposed 
to  stretched  state. 

When  polymers  interfere  with  turbulence,  it  is  because  they  distort  in 
the  fluctuating  velocity  field  and  eventually  react  back  on  the  flow.  The 
approximately  1200 A  size  of  a  typical  polymer  is  much  smaller  than  the 
Kolmogorov  inner  scale  for  most  turbulent  flows.  Thus  the  polymers  see  an 
effectively  uniform  random  straining  field  from  all  the  eddies  in  the  cascade. 
The  Zimm  time  determines  which  eddies  will  have  a  strong  effect  on  the 
polymers  since  they  will  quickly  reassume  their  equilibrium  shape  in  response 
to  strains  induced  by  an  eddy  of  size  /  unless  the  eddy  turnover  time  is  shorter 
than  Tg.  The  strain  rate  7(/)  (i.e.,  the  inverse  eddy  turnover  time)  for  an  eddy 
of  size  /  within  the  Kolmogorov  theory  is  well-known  to  be 

7(/)  =  eV(*.  (1-8) 

The  only  eddies  in  the  turbulent  cascade  which  will  be  able  to  distort  the 
polymer  shape  are  those  for  which  [7]  [8] 

7(0^*  >  (1  -  V 

i.e.,  those  eddies  with  /  <  /*  where 

/•  =  (cr3)i  ^  (1  -  10) 


It  is  appropriate  to  intfject  here  that  the  observation  that  a  critical 
wall  Reynolds  stress  must  be  exceeded  for  drag  reduction  to  result  is  closely 
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related  to  the  time  criterion  expressed  by  Equation  (1-9).  This  is  easily  seen 
by  noting  that  for  wall  flows,  the  maximal  eddy  shearing  rate  is 
7  =  UHu  =  U»/xv,  where  x„  =  vjU^  is  the  width  of  the  viscous  sublayer. 
Thus,  drag  reduction  will  occur  if  £/,/x„  >  t“^,  which  is  equivalent  to  > 
vjTz.  This  is  equivalent  to  the  condition  that 

T„  =  (,ui  >  ksTipe,,  (1  - 11) 

which  defines  a  critictil  wall  Reynolds  stress  varying  inversely  with  i?®  and 
which  is  independent  of  solvent  density,  viscosity  and  pipe  radius. 

The  length  scale  /*  is  independent  of  the  polymer  concentration,  and  is 
of  order  10“^cm  for  a  typical  c  «  10®CTn*/scc®.  Possible  consequences  of  the 
polymer  distortions  induced  by  eddies  smaller  than  this  scale  but  larger  than 
the  Kolmogorov  inner  scale  will  be  discussed  in  Sections  2  and  3.  Here,  we 
simply  note  that  strongly  stretched  polymers  will  begin  to  react  back  on  the 
flow  on  those  scales  1,1^  <  I  <  I*,  when  the  restoring  forces  per  unit  volume 
are  comparable  to  the  Reynold’s  stresses  in  the  turbulent  fluid[8].  For  “ideal” 
polymers,  whose  equilibrium  size  is  given  by  Equation  (1-3),  the  restoring 
force  /  for  a  single  polymer  is  linear  in  the  distorted  polymer  size  Q[6], 

f  =  KQ  =  ksTQIRl.  (1  - 12) 

Here,  K  is  the  Hookean  spring  constant  discussed  above,  and  Q,  (see  Section 
3)  can  be  defined  in  terms  of  a  hydrodynamic  average  of  a  stretching  vector 
describing  the  distortion  of  each  polymer.  Note  that  the  energy  associated 
with  this  force,  E  =  \KQ^,  is  comparable  to  the  thermal  energy  ksT  when 
Q  =  Rg.  This  result  applies  only  for  weak  distortions  of  polymers.  The  force 
will  be  much  larger  when  the  polymer  is  fully  extended,  i.e.,  for  Q  ~  Na. 
When  excluded  volume  interactions  are  taken  into  account,  there  is  also  an 
intermediate  regime  of  stretching,  Rg  «  Q  «  Na,  for  which[6] 

/  «  kBTQ^IRl.  (1  -  13) 

A  complete  theory  of  polymer  drag  reduction  must  determine  how  the  non¬ 
linear  springs  embodied  in  Equation  (1-13)  alter  the  flow  whidi  stretch  them. 
The  construction  of  such  a  theory  is  a  formidable  challenge;  as  we  shall  see 
in  Section  3,  the  hydrodynamics  of  turbulent  flows  in  which  the  polymers 
behave  as  linear,  Hookean  springs  is  already  quite  nontrivial. 
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1.3  Branched  Polymers  and  Polymerized  Membranes 


The  simple  ideas  sketched  above  atilow  us  to  estimate  whether  the  syn¬ 
thesis  of  polymers  with  novel,  nonlinear  architectures  would  be  useful  in 
turbulent  drag  reduction.  It  is  possible  to  synthesize  branched  polymers, 
for  example,  whose  connectivity  is  illustrated  in  Figure  1-2.  There  is  cur¬ 
rently  considerable  interest  in  the  other  novel  architecture  shown  in  Figure 
1-2,  that  of  polymerized  membranes{9].  These  flexible  “tethered  surfaces” 
assume  configurations  in  a  solvent  not  unlike  crumpled  pieces  of  paper. 


Upon  first  consideration,  branched  polymers  and  polymerized  membranes 
appear  to  be  disadvantageous  compared  with  linear  polymers.  The  key  rela¬ 
tion  is  the  scaling  of  the  radius  of  gyration,  with  the  number  of  monomer 
units,  N: 

Rg  «  a'n*';  (1  —  14) 

here  n  is  the  number  of  nodes  in  a  branched  polymer  or  tethered  surface, 
and  a'  is  proportional  to  where  m  is  the  average  number  of  monomers 
per  nodes  (mn  w  N).  The  exponent,  i/,  is  approximately,  0.5  for  branched 
polymers [10],  and  0.4  for  polymerized  surfaces [9]. 


For  fixed  iV,  the  radius  of  gyration  decreases  as  the  number  of  nodes 
increases,  which  makes  the  polymer  less  effective  in  drag  reduction.  One 
wishes  to  increase  the  Zimm  time  as  much  as  possible  so  as  to  ’activate’  the 
polymers  by  eddy-induced  stretching  at  increasingly  larger  onset  scales  given 
by  (see  Equation  (1-10)),  ig  Yet,  the  Zinun  time  is 

proportional  to  il®,  eind  hence  is  shorter  (for  fixed  N)  for  branched  polymers 
or  tethered  surfaces  compared  to  linear  polymers.  Second,  one  wants  to 
permeate  a  significant  fraction  of  the  fluid  with  polymer  even  though  they 
are  a  negligible  component  by  weight.  The  relevant  figure  of  merit  is  the 
mass  density  of  monomer  inside  the  radius  of  gyration. 


Na  \ri) 


N' 


-3if 


(1  - 15) 


where  Mq  is  the  monomer  molecular  weight  and  Na  is  Avogadro’s  number. 
For  a  fixed  number  of  monomers,  AT,  the  cheapest  additives  will  be  those 


10 


Figure  1-2.  Various  chemical  linkages  which  could  be  used  for  turbulent  drag  reduction;  (a) 
conventional  lineeir  polymer;  (b)  Branched  polymer  characterized  by  a  fixed  ratio 
of  branch  points  to  monomer  units;  (c)  A  polymerized  membrane,  with  the 
internal  connectivity  of  a  flat  two  dimensional  surface.  Such  surfaces  can  be  made 
more  flexible  by  connecting  the  nodes  with  polymer  chains,  as  in  a  two 
dimensional  gel. 


which  produce  the  most  tenuous  structures,  i.e.,  those  with  the  lowest  value 
of  p  or  the  largest  v.  Not  surprisingly,  linear  polymers  are  also  most  efficient 
from  this  point  of  view. 

However,  branched  polymers  or  tethered  surfaces  have  the  potential  ad¬ 
vantage  that  they  will  be  less  susceptible  to  degradation  at  high  Reynolds 
numbers.  When  linear  polymers  are  fully  stretched  in  the  extensional  flow 
fields,  they  break  (usually  near  the  center)  and  become  less  effective  in  reduc¬ 
ing  drag.  To  appreciate  this  effect,  let  us  model  the  stretched  linear  polymer 
as  a  set  of  spheres  of  radius  h  connected  by  springs  with  spring  constant  Jfc, 
as  shown  in  Figure  1-3.  We  will  assume  that  the  polymer  is  stretched  lin¬ 
early  along  the  x-direction  by  an  extensional  flow  field,  =  7X,  as  might  be 
created  locally  by  turbulent  eddies,  and  that  the  flow  field  is  centered  on  the 
middle  of  the  polymer.  If  only  the  spheres  interact  hydrodynamically  with 
the  flow  field,  the  net  force  on  the  n**  sphere  is 


Fn  =  -  hSn  -  SwTihyXn,  (1  -  16) 


where  n  =  0  is  the  center  of  the  polymer,  Xn  is  the  position  of  the  sphere, 
—  In  —  ®  is  the  extension  of  the  n*^  spring  from  its  equilibrium 
length  X,  and  p  is  the  viscosity  of  the  fluid.  In  equilibrium,  F  =  0.  One 
immediately  observes  that  >  8n;  in  other  words,  springs  closer  to  the 
center  of  the  polymer  are  stretched  more.  If  we  take  the  continuum  limit, 
repladng  n  by  the  continuous  variable  s.  Equation  (1-16)  becomes 


d8{s) 

ds 


-Ox(5), 


(1  - 17) 


where  a  =  QirijbjJk.  Near  the  center  of  the  spring,  ds  a  dx,  so  Equation  (1- 
16)  has  the  solution  S  {a/2){L^  —  x’),  indicating  a  parabolic  falloff  in  the 
tension  as  one  moves  away  from  the  center  of  the  polymer.  This  argument 
may  explain  why  experiments  on  polymer  degradation  at  ultrahigh  strain 
rates  indicate  breaking  near  the  middle  of  the  chain. 


The  branched  polymer  illustrated  in  Figure  1-3  (or  a  tethered  surface) 
may  alleviate  the  degradation  problem.  Let  us  first  suppose  that  the  branched 
polymer  is  arranged  so  that  the  end  with  fewer  branches  is  furthest  away  from 
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O  o  o  (3XjUULL1XD^JUIJ!JUU0^ 


Figure  1-3.  Polymer  ball-and-spring  model  in  a  high  Reynolds  number  extensional  flow 
field,  V  =  The  center  of  the  polymer  is  at  x  =  0  and  all  springs  are 
stretched  nearby  along  the  x-axis.  Fs  is  the  stoker  drag  on  the  Nth  ball, 

Fs  =  eTijb-yXN-  [(a)  Linear  polymer  (b)  Bremch  polymer.] 
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the  center  of  the  extensional  flow  field  (see  remarks  below).  The  tension  in¬ 
creases  as  one  proceeds  inwards  from  the  end  of  the  polymer  until  one  reaches 
a  node.  At  the  node,  there  are  two  (or  more)  polymer  strands  which  share 
the  load;  that  is,  the  coefficient  of  the  first  term  in  Equation  (1-15)  is  multi¬ 
plied  by  two  (or  more),  so  now  S„-i  may  be  less  than  6„.  Adding  more  nodes 
as  one  proceeds  closer  to  the  center  of  the  flow  field  continues  to  limit  the 
tension,  and,  hence,  avoids  degradation  at  higher  Reynolds  numbers.  Given 
that  the  radius  of  gyration  decreases  as  the  number  of  nodes  increases,  one 
wishes  to  have  the  minimum  number  of  nodes  necessary  to  avoid  degradation 
at  the  design  Reynolds  number. 

In  simple  models,  branched  polymers  begin  from  linear  chains  that  grow 
out  from  a  central  core  and  then  bifurcate  after  some  fixed  average  number 
of  monomer  units;  the  branches  then  bifurcate  after  growing  and  extending 
a  similar  number  of  monomer  units,  etc.  If  such  a  branched  polymer  were 
placed  in  an  extensional  flow  field  with  the  central  node  in  the  middle,  the 
effect  would  be  the  opposite  of  that  described  above.  The  tension  on  the 
monomers  in  the  many  distant  limbs  would  be  focused  on  the  central  node. 
This  may  not  be  a  real  problem  since  this  would  merely  cause  the  original 
branches  to  snap  from  the  central  node,  in  which  case  each  fragment  would 
be  similar  to  the  cartoon  in  Figure  1-3.  However,  it  remains  unclear  whether 
the  fragments  will  align  themselves  in  the  flow  field  so  that  the  branch  that 
was  central  originally  moves  to  the  farthest  part  in  the  flow  field. 

Polymerized  surfaces  may  be  better  candidates  in  that  their  behavior  is 
more  predictable.  We  might  envisage  the  surface  as  a  square  grid  of  polymer 
chains  meeting  at  nodes.  In  an  extensional  flow  field,  the  grid  will  deform 
into  an  extended,  nearly  linear  configuration  with  the  comers,  say,  furthest 
from  the  center  of  the  extensional  flow  field.  As  with  the  branched  polymers, 
the  nodes  in  the  grid  can  limit  the  tension  on  segments  near  the  center  of  the 
flow  field.  The  disadvantage  of  polymerized  surfaces  compared  to  branched 
polymers  is  that  1/  is  smaller.  However,  if  the  grid  is  rather  sparse  (few 
nodes),  this  may  not  be  so  significant. 

In  conclusion,  it  appears  that  the  application  of  polymers  with  novel 
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architectures  is  a  possible  mechanism  for  extending  polymer  drag  reduction 
to  flows  with  very  high  Reynolds  number.  The  addition  of  crosslinks  helps 
the  ploymers  resist  degradation,  but  also  means  that  greater  concentrations 
of  polymers  are  necessary  to  obtain  the  same  drag  reduction. 
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2  THEORIES  OF  DRAG  REDUCTION 


2.1  An  Overview  of  The  Problem 


Drag  reduction  by  polymers  in  turbulent  flows  is  an  extremely  compli¬ 
cated  problem.  It  combines  the  complexity  of  turbulent  flow  (difficult  even 
for  a  Newtonian  fluid)  with  the  problems  of  polymer  physics;  their  combina¬ 
tion  changes  the  character  of  the  turbulence  and  leads  to  a  yet  more  complex 
and  difficult  problem. 

A  naive  attempt  to  understand  the  effects  of  polymers  on  drag  could  go 
astray  in  any  one  of  several  ways.  Dissolved  polymers  increase  the  viscos¬ 
ity  of  polymer  solutions,  particularly  if  the  polymer  molecules  are  extended; 
it  might  therefore  be  (incorrectly!)  inferred  that  they  would  increase  drag 
(albeit  slightly),  just  as  would  an  increase  in  the  viscosity  of  a  pure  New¬ 
tonian  solvent,  even  at  high  Reynolds  numbers.  The  fact  that  in  turbulent 
flows  the  observed  effect  is  in  the  opposite  direction  is  sufficient  to  establish 
the  subtlety  of  the  physics  involved.  So  long  as  the  flow  remains  laminar, 
the  naive  drag  increase  is  expected  to  be  correct;  however  in  laminar  flow 
the  polymers  usually  remain  coiled  so  this  effect  would  be  very  small.  It 
instead  might  be  supposed  that  polymers  could  have  no  effect  at  all  on  drag, 
because  the  rate  at  which  viscous  work  is  done  is  determined,  in  a  stan¬ 
dard  Kolmogorov  cascade  picture,  by  the  largest  scale  velocities  of  the  flow. 
However,  polymer-enhanced  viscous  stresses  are  significant  only  on  the  small 
scales;  which  satisfy  the  Zimm  time  criterion  u,  <  V{t)lt.  Larger  eddys,  for 
which  V{l)lt  <  Ug,  are  unaffected  by  the  presence  of  polymers.  However  it 
must  be  noted  that  momentum  transport  in  bounded  shear  flows  is  critically 
sensitive  to  the  dynamics  of  such  small  scale  eddys  inunediately  outside  the 
viscous  sublayer. 

It  is  evident  from  the  data  (Figure  1-1)  that  polymer  additives  do  not 
affect  the  onset  of  turbulence.  This  is  probably  due  to  the  fact  that  the  un¬ 
stable  sheared  flows  which  commonly  proceed  the  onset  of  turbulence  rotate 
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polymer  chcuns  but  do  not  stretch  them.  Their  effect  is  felt  in  turbulent 
flow,  which,  for  these  data,  occurs  when  Re  >  2500.  At  such  high  Reynolds 
numbers  turbulence  is  expected  to  be  fully  developed,  and  manifest  a  cas¬ 
cade  of  turbulent  energy.  Yet,  as  just  discussed,  the  properties  of  such  a 
cascade  (except  for  the  location  in  k-sp&ce  of  its  inner  scale)  are  independent 
of  the  magnitude  of  the  viscosity.  These  properties  include  the  dissipation 
rate,  which  is  determined  by  the  flow  at  the  outer  (largest)  scale,  and  whose 
volume  average  is  proportional  (by  conservation  of  energy)  to  the  global  flow 
drag  or  friction  factor. 

The  (Kolmogorov)  cascade  model  must  therefore  be  modified  for  polymer 
solutions.  One  possible  modification  is  simply  in  the  numerical  coefficients. 
The  Kolmogorov  analysis  is  entirely  dimensional,  and  cannot  determine  the 
coefficients.  It  is  plausible  that  the  non-Newtonian,  anisotropic,  and  hys- 
teretic  nature  of  the  viscosity  of  a  polymer  solution  changes  the  effective 
values  of  the  coefficients.  This  possibility  can  be  tested  by  detailed  numeri¬ 
cal  simulation  of  turbulence  in  such  fluids. 

A  second  possibility  is  that  the  spatial  variation  in  the  effective  viscos¬ 
ity  of  polymer  solution  pipe  flow  is  essential.  In  this  scenario,  advanced  by 
Lumley(7],  polymers  in  the  boundary  layer  are  thought  to  be  stretched  and 
extended  by  inertial  range  eddys!  As  a  result  of  the  ensuing  coil-stretch  tran¬ 
sition,  the  viscosity  increases  dramatically.  It  is  important  to  keep  in  mind 
that  this  buffer-layer  region  of  enhanced  viscosity  is  located  immediately  next 
to  the  viscous  sublayer,  and  extends  to  a  distance  from  the  wall  x  = 
thus  encompassing  the  eddys  that  satisfy  the  Zimm  criterion.  As  a  result 
of  the  locally  increased  viscosity,  transport  of  stream- wise  momentum  to  the 
wall  decreases  (due  to  small  scale  eddy  damping)  and  drag  is  thus  reduced. 

A  third  possibility  is  that  the  Kolmogorov  cascade  is  truncated,  or  is 
even  inapplicable  because  energy  propagates  up  the  cascade  (to  longer  wave¬ 
lengths)  instead  of  down  it,  violating  Kolmogorov’s  central  assumption. 

An  important  question  not  directly  addressed  by  experiments  to  date  is 
whether  polymer  extension  is  required  for  drag  reduction.  Most,  but  not 
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all,  theories  would  predict  that  this  is  so.  The  prediction  could  be  tested  in 
experiments  with  sufficiently  concentration  solutions  of  short-chain  polymers. 
At  the  onset  of  turbulence  the  strain  rates  will  be  insufficient  for  them  to 
undergo  the  coil-stretch  transition,  and  most  theories  will  predict  little  drag 
reduction.  At  high  Reynolds  numbers  and  higher  turbulent  strain  rates,  there 
should  be  an  onset  (sudden  or  gradual)  of  drag  reduction. 

2.2  Polymers  in  Turbulent  Flow 

The  behavior  of  dissolved  polymers  in  flow  fields  has  been  the  subject  of 
an  extensive  literature.  For  example,  Henyey  and  Rabin[l  1]  have  calculated 
equilibrium  polymer  conformations  in  stationary  flow  fields.  The  problem  of 
dynamic  polymer  relaxation  to  an  equilibrium  state  in  specified  stationary 
flow  fields  is  considerably  more  complex  (Rabin,  Henyey  and  Pathria)[12]; 
Rabin[13].  The  problem  of  polymers  in  turbulent  flow  fields  is  yet  more  diffi¬ 
cult:  the  flow  field  is  varying  spatially  and  temporally,  and  the  quantitative 
nature  of  these  variations  can  be  specified  only  statistically  (and  even  this 
specification  is  both  difficult  and  uncertain). 

There  are  a  number  of  models  and  an  extensive  literature  on  the  subject 
of  dissolved  polymer  response  to  flow  fields.  Two  simple  limiting  models 
are  the  dumbell  model  (for  example,  as  developed  in  a  later  chapter  of  this 
report)  and  the  affine  deformation  model,  in  which  the  polymer  is  passively 
deformed  following  the  strain  tensor  of  the  embedding  fluid.  It  is  frequently 
argued  that  a  great  degree  of  polymer  extension  (necessary  in  order  that  a 
low  polymer  concentration  produce  significant  effects  on  the  flow)  requires 
that  some  component  of  the  strain  rate  tensor  exceed  the  Zimm  relaxation 
rate  of  the  polymer,  proportional  to  N~^'^  in  a  good  solvent,  where  N  is 
the  number  of  monomers  in  the  polymer  chain.  This  argument  has  been 
challenged  by  RHP[12]  who  present  a  free  energy  argument  which  predicts 
a  dependence  proportional  to  AT"'  *  in  good  solvents.  They  assert  that  this 
latter  form  is  in  better  agreement  with  the  data.  Both  these  arguments  are 
essentially  kinetic,  in  that  they  depend  on  overcoming  the  entropic  barrier 
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to  full  polymer  extension;  in  any  significant  stationary  extensional  flow  field 
the  equilibrium  state  of  a  long  polymer  is  nearly  fully  extended.  This  follows 
from  the  Kramers  potential  describing  the  tendency  of  the  viscous  drag  force 
to  distend  the  polymer,  which  varies  as  a  steeper  dependence  than  the 
entropic  forces  which  tend  to  keep  it  coiled.  For  small  extensions  (coiled 
polymers)  R  is  small  and  the  entropic  forces  dominate.  The  arguments  differ 
in  that  the  RHP  argument  demands  that  there  be  no  entropic  barrier,  while 
the  Zimm  time  argument  demands  only  that  the  entropic  barrier  be  ineffec¬ 
tive  (according  to  an  approximate  argument  based  on  relaxation  rates).  One 
might  question  whether  there  can  be  anything  like  this  sharp  phase  transition 
to  an  elongated  state  when  the  flow  is  highly  turbulent;  phase  transitions  are 
often  smeared  out  in  systems  subjected  to  chaotic,  time-dependent  forces[6]. 


The  problem  of  polymer  extension  in  turbulent  flow  fields  is  particularly 
difficult  and  subtle.  The  principal  axes  of  extension  are  changing  with  a 
characteristic  decorrelation  time  r^.  The  flow  has  a  fluctuating  vorticity  of 
characteristic  magnitude  w,  which  rotates  the  polymers,  and  tends  to  inter¬ 
fere  with  their  extension  (it  is  for  this  reason  that  relatively  little  extension  is 
obtained  in  pure  shear  flow,  for  which  the  magnitude  of  vorticity  equals  the 
largest  component  of  the  rate  of  strain  tensor).  Of  course,  for  Kolmogorov 
turbulence,  w  ~  rj*.  At  high  strain  rates  a  polymer  molecule  will  have  neg¬ 
ligible  entropic  restoring  forces,  and  will  passively  (affinely)  follow  the  flow. 
In  a  stationary  purely  extensional  flow  field  with  the  rate  of  strain  dv/dr  its 
length  will  increase 

£  w  /<,exp  •  (2-1) 


In  the  turbulent  flow  field, 
be 


after  a  single  decorrelation  time  the  extension  will 


£  into  exp 


(2-2) 


where  the  effects  of  vorticity  are  qualitatively  included  by  defining  in  axes 
and  along  a  path  which  follow  the  polymer  molecule. 


It  is  evident  that  the  polymer  extension  depends  very  sensitively  on 


(2-3) 
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As  noted  above,  this  quantity  may  be  estimated,  on  dimensional  grounds,  to 
be  of  order  unity  for  a  turbulent  flow.  It  is  clear  that  its  actual  quantitative 
value  (or  rather,  that  of  the  analogous  quantity  in  a  more  quantitative  the¬ 
ory)  is  of  great  importance,  and  is  a  sensitive  measure  of  the  nature  of  the 
turbulent  flow.  However,  it  does  appear  unlikely  that  significant  extension  is 
possible  in  highly  turbulent  flows. 

2.3  Drag  Reduction 

Very  little  of  the  extensive  literature  on  polymers  in  turbulent  flow  actu¬ 
ally  attempts  to  explain  why  there  should  be  a  reduction  in  drag.  The  classic 
work  of  Lumley[7]  argues  that  if  polymer  molecules  are  extended  by  turbu¬ 
lence  in  pipe  flow,  they  lead  to  a  flow  in  which  two  definitions  of  the  laminar 
boundary  thickness  are  possible  -  one  based  on  the  viscosity  with  unextended 
polymers,  which  should  be  applicable  in  the  laminar  viscous  sublayer,  and 
a  larger  value  calculated  from  the  viscosity  in  the  turbulent  region  with  ex¬ 
tended  polymers.  He  argues  that  between  these  two  distances  from  the  wall 
there  should  be  a  buffer  layer  into  which  the  turbulence  has  difficulty  pene¬ 
trating,  and  thus  through  which  momentum  transfer  is  slow,  thus  reducing 
the  wall  drag.  This  buffer  layer  corresponds  to  the  region  of  enhanced  slip 
(i.e.  good  streamwise  momentum  confinement).  This  model  is  qualitatively 
reasonable,  but  is  not  quantitative  enough  to  test. 

Ryskin[14]  describes  a  “yo-yo”  model  for  affine  polymer  deformation  in 
extensional  flow.  It  appears  to  satisfactorily  describe  convergent  flow  in  a 
cone.  When  applied  to  pipe  flow,  this  model  resembles  Lumley’s;  its  treat¬ 
ment  of  polymer  dynamics  in  a  specified  flow  field  is  quantitative,  but  its 
application  to  turbulence  and  its  treatment  of  drag  reduction  are  partly  phe¬ 
nomenology,  and  partly  semi-empirical  fit  to  the  data. 

Tabor  and  de  Gennes[15]  (see  also  de  Gennes[8]  for  more  detail)  present 
a  very  different  model.  They  assume  affine  polymer  deformation.  In  order 
to  calculate  the  degree  of  extension,  they  assume  that  the  fluid  strain  is  a 
unique  function  of  the  length  scale  of  the  turbulence,  in  analogy  to  laminar 
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cone  flow  in  n  dimensions,  where  n  is  a  non-integral  number  between  1  and 
2.  Each  polymer  molecule  is  assumed  to  be  affected  only  be  eddies  of  a 
single  size.  This  description  is  peculizu’  and  of  dubious  applicability  because 
cone  flow  is  strongly  convergent  and  has  a  large  (negative)  divergence  at 
its  apex,  while  incompressible  pipe  flow  and  turbulence  are  divergence-free. 
Also,  a  given  molecule  is  simultaneously  subjected  to  turbulent  eddies  of 
all  sizes,  as  pointed  out  by  Ryskin[14].  Having  made  these  assumptions, 
they  define  a  new  inner  scale  I**  for  the  turbulent  flow  field,  by  equating 
the  polymer  elastic  energy  density  to  the  Reynolds  stress  of  the  turbulence. 
By  analogy  with  the  early  work  of  Lumley,  but  without  further  quantitative 
detail,  they  argue  that  the  existence  of  this  larger  inner  scale  for  the  polymer 
-  hydrodynamic  turbulent  flow  leads  to  drag  reduction  by  truncation  of  the 
Kolmogorov  cascade. 

The  presence  of  extended  polymers  in  a  turbulent  flow  will  increase  the 
viscous  dissipation  in  eddies  of  a  given  size,  and  thus  increase  the  inner  scale 
size  of  the  flow.  It  is  possible  to  construct  a  model  analogous  to  that  of 
Tabor  and  de  Gennes  in  which  the  new  inner  scale  is  set  by  equating  the 
viscous  dissipation  rate  in  flow  about  the  polymer  molecules  to  the  rate  of 
dissipation  in  the  turbulence: 

=  (2-4) 

where  Cp  is  the  number  density  of  polymers,  R  their  radius  of  gyration,  7 
the  shear  rate,  rj  the  solvent  viscosity  and  p  its  density,  and  U{r)  is  the 
characteristic  velocity  of  eddies  of  size  r.  Following  Tabor  and  de  Gennes  we 
write 

R  =  ARo  (2-5) 

A  =  i^)"  (2-6) 

where  i*  is  the  length  scale  at  which  7(/*)  equals  the  Zinun  relaxation  rate 
(the  criterion  for  afSne  polymer  deformation).  Assuming  Kolmogorov  turbu¬ 
lence  with  dissipation  rate  c  and  estimating  7(r)  =  U{r)/r,  a  little  algebra 
leads  to 


3  A  MODEL  OF  TURBULENT  POLYMER 
HYDRODYNAMICS 

3.1  The  Dumbell  Model 


We  focus  our  attention  on  the  two  end  points  of  a  long  polymer  chain,  and 
replace  its  many  internal  degrees  of  freedom  by  a  stretchable  dumbell  which 
spans  the  points-see  Figure  3-1.  Although  clearly  a  gross  oversimplification, 
this  crude  model  has  proven  useful  in  studies  of  laminar  polymeric  flow[16], 
and  leads  to  nontrivial  modifications  of  homogeneous  isotropic  turbulence. 
The  equations  of  motion  of  the  two  endpoints  are  assumed  to  be 


where  vi  and  V2  are  the  fluid  velocities  at  the  endpoints,  7  w  rj^Rg  is  the 
Stoke’s  drag  coeflScient  for  the  part  of  the  pol”~ier  represented  by  the  end¬ 
point,  and  the  spring  potential  is 

=  (3-2) 

with  K  =  ksT/Rg.  Although  nonlinear  spring  potentials  like  that  repre¬ 
sented  by  Equation  (1-13)  are  more  realistic  in  the  highly  stretched  regime, 
we  shall  for  simplicity  confine  our  attention  to  the  more  tractable  linear  case. 
In  setting  7  ^  Jj,Rg,  we  make  the  usual  assumption[4]  that  the  fluid  in  the 
vicinity  of  the  polymer  chain  is  dragged  along  by  the  chain  motion,  so  that 
the  drag  is  of  order  the  Stokes  force  on  an  object  cf  size  Rg.  The  functions 
Ci(t)  and  ^2(0  thermal  noise  sources,  with  autocorrelation  functions 

<  C.(0Ci(O  >=  2^k,T6ii6{t  - (3  -  3) 

which  serve  to  equilibrate  the  dumbell  when  there  is  no  macroscopic  fluid 
motion. 
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Figure  3-1.  Long  chain  polymer  with  endpoints  Fj  and  r2  modeled  as  a  stretchable  dumbbell. 
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We  imagine  that  the  fluid  is  filled  with  a  dilute  solution  of  these  dumbells 
and  pass  to  center-of-mass  amd  relative  coordinates  for  the  two  ends, 


R  =  5(f.  +  fi) 

q  =  ri-  rj, 

in  terms  of  which  the  Langevin  Equation  (3-1)  may  be  written 

f  =  v+^tC.W+««)l 
f  =  +  i(C.W-W01. 


(3-4) 


(3-5) 


Here  v  is  the  external  velocity  field  evaluated  at  the  center-of-mass  position, 
and  we  have  expanded  vi  —  in  gradients  of  the  velocity  field  at  this  point 
to  get  Equation  (3-5).  Standard  methods[l7]  then  lead  to  a  Fokker-Planck 
equation  for  the  distribution  function  /(fi,g,<)  describing  the  midpoints  and 
end-to-end  distances  of  the  ensemble  of  dumbells,  namely 


(3-6) 


%■  -  s 
- 1 


where  /  is  normalized  so  that 


(3-7) 


The  diffusion  constants  Do  and  D,  arise  from  the  Langevin  noise  sources, 
and  are  given  by 


Do  =  kBTj2'i 
D,  =  2kBTh. 


(3-8) 


In  a  steady  homogeneous  shear  flow  parametized  by  the  constant  matrix 


diVj  — 


(3-9) 


one  can  check  that 


=  exp[-{u{^/kBT)  +  aijqiqj] 


(3  -  10) 


25 


is  a  steady  state  solution  of  (3-6). 


Hydrodynamic  equations  follow  from  taking  various  moments  of  Equation 
(3-6).  The  concentration  c(r,  t)  of  dumbells 

c(r,i)  =  /A/(n9,<)  (3-11) 

for  example,  obeys  the  equation 

dtc  +  {v  -  V)c  =  V=c.  (3  -  12) 

An  initially  inhomogeneous  dumbell  distribution  subjected  to  homogeneous 
isotropic  turbulence  will  be  mixed  to  uniformity  in  a  few  eddy  turnover  times. 
Henceforth,  we  neglect  the  complication  of  inhomogeneities  in  the  polymer 
concentration,  and  simply  set 

c(r,  t)  =cp  const.  (3  —  13) 

It  is  also  illuminating  to  take  the  first  moment  of  Equation  (3-6)  and  consider 
the  dynamics  of  the  average  end-to-end  distance 

Q(.?,t)  =  J  (3-14) 

which  satisfies 

dt0  +  {v  .  V)(5  =  ((?  •  V)v  -  2a;,(5  -|-  P,  V*(5  (3  -  15) 

with  u)x  =  K/'y.  In  addition  to  the  usual  streaming  and  diffusive  terms, 
there  is  a  line  stretching  term  (0  •  V)t>.  This  tendency  of  velocity  gradients 
to  stretch  the  dumbells  is  resisted  by  the  term  — 2u;,Q,  where  u;,  is  of  order 
the  reciprocal  of  the  Zimm  relaxation  time  discussed  in  Section  1.0. 

Equation  (3-15)  is  in  fact  inappropriate  for  the  problem  at  hand;  to  see 
this,  note  first  that,  as  shown  in  Figure  3-2,  many  different  polymers  will 
interact  with  a  given  volume  of  fluid.  Any  attempt  to  define  a  hydrodynamic 
field  0(rjt)  by  averaging  polymer  stretching  vectors  over  this  volume  will 
depend  on  purely  arbitrary  conventions  such  as  whether  we  set  q  =  ri  — 
or  q  =  fi  —  ri  for  a  given  dumbell.  If  one  eliminates  the  ambiguity  by,  say, 
using  polymers  which  are  deuterated  at  one  end,  the  configuration  q  will  be 
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Figure  3-2.  Elongated  dumbbells  in  a  hydrodynamic  averaging  volume,  together  with  the 
microscopic  stretching  vectors  {q®}.  The  appropriate  order  parameter 
describing  the  state  of  alignment  is  =  <q“iq“j>.  The  sense  of  the  arrows 
depends  on  a  arbitrary  convention,  which  drops  out  with  this  definition  of 
the  order  parameter. 
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as  likely  as  —q,  and  the  hydrodynamic  average  Q{f,  t)  will  be  zero,  even  if 
the  stretched  polymers  are  aligned!  The  same  problem  appears  in  nematic 
liquid  crystals,  which  are  fluids  of  aligned  rod-shaped  molecules  with  inver¬ 
sion  symmetry[18].  As  in  the  case  of  nematics  the  appropriate  hydrodynamic 
field  is  a  tensor. 

*)  =  J  g,  0»  (3-16) 

which  is  manifestly  invariant  under  q  —*  —q.  The  equation  of  motion  which 
follows  from  Equation  (3-6)  is 

d,Q,i  +  {S-  ^)Q,J  =  +  D.V^Qii  +  — <0-. 

\z  - 17) 

In  addition  to  streaming,  stretching  and  diffusive  terms  similar  to  those  ap¬ 
pearing  in  Equation  (3-15),  there  is  a  thermal  forcing  contribution  {‘ikBTf'y)Sij. 
This  term  balances  against  the  Zimm  relaxation  term  — in  a  uniform 
quiescent  fluid  to  give  the  expected  equilibrium  solution  of  (3-17),  namely 

=  =  (3-18) 

corresponding  to  an  isotropic  distribution  of  polymers  with  size  R,. 

The  tensor  Qij(r,t)  determines  how  the  stretched  polymer  distribution 
reacts  back  on  the  fluid.  Taking  over  the  analysis  for  laminar  polymeric  flow 
in  Reference  [16],  we  find  that  the  equation  of  motion  for  the  (incompressible, 
V  •  u  =  0)  velocity  field  is 

’  dv’ 

P  -  ~^iP  +  (3  -  19) 

The  new  term,  cpKdjQij,  is  proportional  to  the  concentration  of  polymers 
and  represents  a  contribution  to  the  fluid  stress  tensor  of  the  form 

S<Tii  =  CpKQij.  (3-20) 

As  explained  in  [16],  it  arises  from  the  forces  exerted  on  the  fluid  by  dumbells 
which  pierce  the  walls  of  the  hydrodynamic  averaging  volume. 

Equations  (3-17)  and  (3-19)  comprise  a  nontrivial  hydrodynamical  model 
of  polymeric  fluids  which  would  be  very  interesting  to  study  at  high  Reynold’s 
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numbers  via  computer  simulations.  Except  for  its  tensorial  character  (Qij  has 
six  independent  components),  Equation  (3-17)  is  like  the  equation  which  de¬ 
scribes  the  stretching  of  the  magnetic  field  in  magnetohydrodyncimics  (MHD)(19]. 
We  expect  that  numerical  codes  developed  for  MHD  could  be  readily  adapted 
to  study  polymer  hydrodynamics  (PHD!).  Although  nonlinear  dumbell  springs 
are  easily  incorporated  into  the  velocity  equation,  such  nonlinearities  lead  to 
an  infinite  hierarchy  of  equations  for  moments  involving  g,  which  are  less 
straightforward  to  simulate.  When  Qij  oc  Sij,  as  in  Equation  (3-18),  the 
polymeric  effects  can  be  incorporated  into  an  additional  (osmotic)  contribu¬ 
tion  to  the  total  pressure.  More  generally,  however,  we  must  decompose  Qij 
into  diagonal  and  traceless  parts, 

Qij{r,  t)  =  Qo{?,  t)Sij  +  Qij{r,  t)  (3-21) 

where  T rQij  =  0,  and  observe  that  Qij  contributes  in  a  nontrivial  way  to  the 
dynamics  of  an  incompressible  fluid.  The  isotropic  and  traceless  parts  of  Qij 
are  coupled  together  in  Equation  (3-17)  by  the  stretching  terms. 


Some  insight  into  the  behavior  of  turbulent  polymeric  fluids  follows  from 
the  observation  that,  in  the  absence  of  dissipation  and  thermal  forcing,  only 
the  sum  of  the  fluid  kinetic  energy  and  polymer  potential  energy  is  conserved. 
Indeed,  it  is  str^ughtforward  to  show  from  Equations  (3-17)  and  (3-19)  that 
if 

Etot  —  2  J (3  -  22) 


then 


dEjoT 

~~dir 


=  -J  <^r[pu{diVjf  +  4u;,Q*jfc]  (3  -  23) 


where  is  the  volume  of  the  sample.  We  imagine  that  the  m^lcroscopic 
forcing  of  the  kinetic  degrees  of  freedom  has  been  turned  off  and  ask  how  the 
energy  decays.  Under  circumstances  such  that  the  polymers  are  only  slightly 
stretched,^  cpKQkk  is  a  negligible  fraction  of  the  kinetic  energy  and  the 
conventional  Kolmogorov  cascade  picture  should  be  at  least  approximately 
correct.  We  argued,  however,  in  Section  1.0  that  eddies  with  characteristic 
size  i  <  £*  =  (c/w^)!  would  produce  effectively  uniform  random  shears 
which  would  strongly  distort  the  polymers.  The  same  criterion  emerges  from 
comparing  the  Kolmogorov  estimate  of  the  stretching  matrix  (diVj)  with 
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the  Zimm  rate  — w,  on  the  right-hand  side  of  Equation  (3-17).  Because  of 
their  coiled  equilibrium  state,  polymers  are  capable  of  enormous  extensions 
and  their  potential  energy  can  become  a  significant  fraction  of  the  total 
energy  in  the  flow.  One  “viscoelastic”  scenario[8],[15]  is  that  kinetic  and 
potential  energy  are  exchanged  between  the  fluid  and  polymer  degrees  of 
freedom  over  a  range  of  scales  td  <  I  <  I**  <  I*  in  a  way  which  effectively 
truncates  the  propagation  of  the  kinetic  energy  down  to  the  Kolmogorov 
inner  scale.  Energy  is  then  destroyed  at  a  slower  rate  by  the  dissipative 
terms  on  the  right-hand  side  of  Equation  (3-23),  because  the  average  value 
of  pvipiVjY  decreases.  If  this  effect  of  adding  polymers  is  larger  than  the 
increased  dissipation  caused  by  the  Zimm  term  lOzQkkt  energy  will  decay  at 
a  slower  rate  and  the  drag  will  go  down. 

Tests  and  elaborations  of  these  ideas  via  spectral  closure  methods  and 
direct  numerical  simulations  would  clearly  be  very  desirable.  The  scenario  of 
Tabor  and  DeGennes  described  in  Section  2.0,  in  particular,  could  be  tested 
explicitly.  The  “viscoelastic”  range  of  length  scales  they  propose  bears  some 
sinnularity  to  a  range  of  energy  equipartition  between  magnetic  and  kinetic 
degrees  of  freedom  which  occurs  in  models  of  turbulent  MHD,  mediated  by 
Alfven  waves[20].  Elquations  (3-17)  and  (3-19)  do,  in  fact,  admit  “polymer 
wave”  solutions  which  are  similar  to  Alfven  waves.  To  see  this,  we  first 
imagine  that  random  stretching  by  the  turbulent  velocity  field  on  scales  i  < 
has  momentarily  aligned  the  dumbells,  as  in  Figure  1-3.  This  is  like  the 
generation  of  a  local  region  of  intense  magnetic  field  in  MHD.  We  then  look 
for  excitations  out  of  this  state  of  alignment  which  couple  to  the  velocity 
field.  We  write  the  tensor  Qij  as 

Qij  «  iQoi  +  SQi){Qoi  +  6Qi)  (3-24) 

^  QoiQoj  d"  Qoi^Qj  d"  ^QiQoj 

with  Qo  ±  6(5*,  and  look  for  normal  modes  of  (3-17)  and  (3-19)  of  the  form 

(3-25) 

if  = 

with  6Q  •  k  ^  Vo  '  k  =:  0.  For  wavevectors  k  such  that  /J*  <  it  <  (^rf)“*,  we 
can  neglect  both  the  viscous  damping  and  the  relaxation  terms  in  Equation 
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(3-17)  and  find  the  characteristic  frequency 


w(k)  =  ./5e^(4  ■  fy.  (3  -  26) 

V  Po 

We  have  assumed  that  the  local  alignment  produces  a  Qo  which  greatly 
exceeds  Rg,  so  that  the  forcing  term  in  Equation  (3-17)  can  also  be  neglected. 

For  the  above  calculation  to  be  meaningful,  the  frequency  Equation  (3-26) 
must  be  higher  that  the  inverse  turnover  time  of  the  eddy  which  produced 
the  region  of  alignment,  i.e.,  we  require 

w{k)>cika.  (3  —  27) 

In  the  region  of  wavevectors  where  this  happens,  kinetic  energy  will  convert 
into  potential  energy  during  the  lifetime  of  the  aligned  region,  and  the  con¬ 
ventional  Kolmogorov  cascade  will  be  modified  or  terminated.  To  calculate 
where  this  region  occurs,  we  must  evaluate  Ekjuation  (3-26),  and  hence  es¬ 
timate  the  root  mean  square  stretching  Qo  on  a  scale  £  ~  fc"*  due  to  all 
relevant  larger  scales.  Following  heuristic  methods  developed  for  MHD  tur- 
bulence[20],  we  write  this  as 

Ql{k)  =  ^l^^\kk‘^u(k)  (3-28) 

where 

QiiW  =  /  (3  -  29) 

and  where  the  lower  limit  arises  because  scales  t  >  t,  have  no  effect  on  the 
polymers.  Knowing  the  function  Qoi^)  is  equivalent  to  knowing  the  exponent 
n  in  Elquation  (2-8).  Simulations  of  Elquations  (3-17)  and  (3-19)  would  help 
test  this  picture,  and  determine  the  exponent  n. 

3.2  Polymer  Hydrodynamic  Turbulence 

Having  described  the  dumbell  modd  of  polymer  hydrodynamics,  we  now 
turn  to  a  discussion  of  the  dynamics  of  turbulent  polymer  solutions.  While 
bounded  pipe  and  channel  flows  are  of  ultimate  interest,  we  here  focus  on 
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the  simpler  (yet  quite  formidable)  problem  of  understanding  the  physics  of 
homogeneous  polymer  hydrodynamic  turbulence  produced  by  the  flow  of  the 
solution  thru  a  fine  grid  (hereafter  referred  to  as  grid  turbulence).  Our  ap¬ 
proach  is  to  exploit  the  formal  analogy  between  polymer  hydrodynamics  and 
magnetohydrodynamic  turbulence,  the  essential  physics  of  which  is  sununa- 
rized  in  Appendix  A.  As  noted  previously,  the  polymer  hydrodynamic  ex¬ 
tension  field  Q  is  analogous  to  the  magnetic  field  B,  so  that  polymer  elastic 
waves,  which  propagate  at  =  {CpK  are  the  formal  counter¬ 
part  of  Alfven  waves,  which  propagate  at  Va  =  An  important 

difference  between  the  two  systems  occurs  in  the  scale  dependence  of  the 
dissipation.  In  the  case  of  magnetohydrodynamics,  magnetic  energy  is  dissi¬ 
pated  by  resistive  diffusion,  so  that  characterizes  the  damping  rate  of  B*. 
In  contrast,  dissipation  of  elastic  energy  occurs  at  the  Zimm  relaxation  rate 
w,,  which  is  scale  independent.  Hence,  while  the  magnetic  Reynolds  number 
^A/(^)  =  fy{f)h  decreases  with  eddy  scale  size  £,  the  analogous  ’elastic’ 
Reynolds  number  Re(^)  =  V(^)/A»;,  actually  increases  as  £  decreases,  until 
the  rather  feeble  effects  of  polymer-concentration  diffusivity  assert  themselves 
at  microscopic  scales  (i.e.  typically  <  v,  where  u  is  the  fluid  kinematic 
viscosity).  Hence,  dissipative  effects  associated  with  polymer  dynamics  are 
most  important  for  scales  which  marginally  satisfy  the  Zimm  criterion  (i.e. 
scales  i  such  that  V{£)/£  *^0;,),  and  decrease  in  importance  throughout  the 
remainder  of  the  inertial  range. 

In  the  production  of  grid  turbulence,  solution  flows  thru  a  hot  wire  grid 
of  mesh  spacing  L  at  velocity  V^.  Turbulence  then  evolves  (rapidly)  as  the 
fluid  flows  downstream  of  the  grid  (i.e.  see  Figure  3-3).  The  downstream 
evolution  of  pol}rmer  hydrodynamic  turbulence  may  be  considered  to  occur 
in  two  sequential  stages: 

1.  the  elastization  phase[21],[22],[23]  —  a  transient  phase  during  which  the 
heretofore  hydrodynamically  ‘passive’  polymer  coils  are  stretched,  so 
that  some  of  the  mechanical  energy  of  inertial  range  eddys  is  converted 
into  polymer  elastic  energy. 
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Turbulence 


Vo 


Figure  3-3.  The  Production  and  Evolution  of  Grid  Turbulence. 
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2.  the  viscoelastic  turbulence  phase{20],{24]  —  the  stage  of  fully  developed 
viscoelastic  turbulence  during  which  the  coupling  of  hydrodynamic  ed- 
dys  to  elastic  waves  modifies  the  spectrum  and  energy  dissipation  rate 
of  the  turbulence.  Throughout  this  discussion,  it  is  important  to  keep 
in  mind  that  the  turbulence  is  forced  at  the  constant  rate  c  =  V^/L. 
Hence,  e  should  be  thought  of  as  a  forcing  rate,  rather  than  as  a  “dis¬ 
sipation  rate,”  its  usual  meaning.  The  actual  dissipation  rate  will  be 
computed  using 

=  -  /  <Pr  £)»  +  2u,,  0^} .  (3  -  30) 

Also,  it  is  understood  that  the  turbulence  dynamics  are  modified  only 
on  scales  for  which  the  Zimm  criterion  is  satisfied  (i.e.  for  which  V{£)/£  > 
w,). 


In  the  elastization  phase,  passive,  initially  coiled  polymers  are  stretched 
by  turbulent  eddys.  Since 


^  =  2-2K-u-.£ 


(3-31) 


it  follows  that  the  rate  of  polymer  stretching  due  to  eddys  of  scale  £  is  given 

(^)  (3-32) 

Stretch  * 

Thus,  perturbations  in  the  elasticity  field  Q  grow  exponentially  (on  all  scales 
for  which  the  Zimm  criterion  is  satisfied)  at  the  rate  V{£)/£  ~  e  For 
Kolmogorov  turbulence,  small  scale  elasticity  field  perturbations  grow  most 
rapidly.  Of  course,  for  £<  £4  the  growth  rate  drops,  as  such  scales  fall  into 
the  hydrodynamic  dissipation  range.  This  process  of  eddy-stretching  induced 
growth  of  elasticity  perturbations  continues  until  the  back-reaction  of  elastic 
energy,  induced  by  tension  in  the  polymer  strands,  becomes  comparable  to 
the  mechanical  energy  of  the  stretching  eddy.  Thus,  elastic  energy  growth  is 
limited  by  equipartition,  so  that 


^rms  —  PoV^{£xl 


(3  -  33) 


Equation  (3-33)  assumes  £z  <  L,  where  ~  u;,,  i.e.  the  largest  acti¬ 

vated  scale  of  Q. 

In  the  ensuing  phase  of  fully  developed  viscoelastic  turbulence,  two  types 
of  excitations  co-exist  in  the  polymer  solution.  The  first  are  (low  frequency) 
fluid  eddys,  with  w(^)  V{t)lt.  The  second  are  (high  frequency)  elastic 
waves,  with  u){t)  Here  Ve  is  the  shear  elastic  wave  velocity  computed 

using  (?rms,  i-e. 

V|  =  C,K  Q\^h..  (3  -  34) 

As  Ve  contains  contributions  from  all  eddys  with  £  <  tz^  Ve  >  V{t)  for 
£  <  £z.  This  state  of  viscoelastic,  polymer  hydrodynamic  turbulence  is  closely 
analogous  to  magnetohydrodynamic  turbulence,  the  basic  properties  of  which 
are  summarized  in  Appendix  A.  As  in  the  case  of  MHD  turbulence,  elastic 
waves  tend  to  impede  inertial  range  transfer,  since  eddy  interaction  can  occur 
only  when  two  neighboring  elastic  waves  propagating  in  opposite  directions 
collide[14],  generating  a  low  frequency  virtual  mode  which  interacts  with  the 
fluid  eddys  (see  Figure  3-4).  Consequently,  as  in  MHD,  the  inertial  range 
transfer  rate  for  viscoelastic  turbulence  is  reduced[20],  i.e. 

Hence,  self-similarity  of  energy  transfer  implies  V{£)  —  and 

E{k)  =  {tVEfl'^k-^/\  (3-36) 

where  Ek  Eq  ^(^)*  Note  that  the  inertial  range  spectrum  of  PHD 
turbulence  is  A:"®/*  (as  in  MHD)  rather  than  the  familiar  k~^l^  of  the 
Kolmogorov  theory.  Furthermore,  the  dissipation  scale  of  PHD  is  also  dif¬ 
ferent  from  that  of  Kolmogorov  turbulence  and,  following  the  arguments  of 
Appendix  A,  is  easily  seen  to  be 

li,  =  (3  -  37) 

Perhaps  the  most  basic  question  one  can  ask  about  PHD  grid  turbu¬ 
lence  is:  does  adding  the  polymers  increase  or  decease  the  rate  of  energy 
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Figure  3-4.  Eddy  —  Elastic  Wave  Interaction. 
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dissipation?!  This  question  can  now  be  answered  by  simply  computing  the 
difference  of  the  energy  dissipation  rates  for  PHD  turbulence  and  ordinary 
hydrodynamic  turbulence,  i.e.  by  calculating: 

Ac  =  {yk^Ekik)  +  2ij,EQ{k))  (3-38) 

-  uPE{k). 

Here,  Ek  and  Eq  are  the  (equipartitioned)  mechanical  and  elastic  energies 
which  add  to  the  result  given  in  Equation  (3-36),  kd^  =  idpjkd  = 

(the  Kolmogorov  microscale),  k^  =  and  E{k)  =  Note  that  the 

lower  bound  of  the  range  ot  integration  is  k^  in  each  case,  since  only  those 
scales  on  which  the  polymers  are  ’activated’  can  exhibit  altered  dynamics 
and  dissipation  rates.  Straightforward  arithmetic  indicates  that: 

Ac>  0  (3  -  39) 

so  that  the  net  dissipation  rate  is  increased!  The  key  effect  in  determining 
that  Ac  >  0  is  that  Zimm  relaxation- mode-induced  dissipation  is  not  ‘small’ 
for  scales  with  £  «  (i.e.  see  Figure  3-5),  as  is  viscous  dissipation.  Alter¬ 

natively  put,  energy  must  be  dissipated  in  fluid-polymer  friction  drag  during 
the  process  of  stretching  the  polymer.  It  should  also  be  mentioned  that  more 
detailed  polymer  models  would  consist  of  a  spectrum  of  Zimm  relaxation 
modes,  each  with  a  corresponding  relaxation  frequency.  In  that  case,  one 
might  expect  increased  Zimm  dissipation  losses,  since  a  broad  spectrum  of 
dissipation  modes  (rather  than  just  one)  would  resonate  with  the  fluid  eddys. 
If  a  sufficient  number  of  relaxation  modes  were  excited  the  cascade  could  be 
effectively  truncated.  Finally,  we  note  that  the  prediction  that  Ac  >  0  should 
be  amenable  to  experimental  examination. 

At  first  glance,  the  result  of  Equation  (3-39)  seems  disturbing,  since  drag 
is  predicted  to  increase!  However,  in  interpreting  the  result  one  must  keep  in 
mind  the  important  differences  between  grid  turbulence  and  pipe  flow  turbu¬ 
lence.  In  particular,  the  results  of  the  above  analysis  suggest  that  in  a  pipe 
flow  a  viscoelastic  buffer  layer  will  form  in  the  region  xy  <  x  <  x,  =  f/./w, 
of  the  turbulent  boundary  layer  (see  Figure  3-6).  Here  xy  =  is  the 
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Figure  3-6.  The  Viscoelastic  Boundary  Layer, 
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width  of  the  viscous  sublayer.  Within  this  viscoelastic  buffer  layer,  cill  ed- 
dys  satisfy  the  Ziirun  activation  criterion,  and  thus  one  can  expect  a  mod¬ 
ified  inertial  range  spectrum  and,  most  importantly,  enhcinced  dissipation 
in  the  layer.  The  enhanced  dissipation  in  turn  results  in  reduced  momen¬ 
tum  transport,  and  thus  in  enhanced  slip  and  drag  reduction.  We  note  that 
this  scenario  is,  in  effect,  a  variation  on  the  ‘Newtonian  plug  and  buffer 
layer’  scen2irio  of  Lumley,  with  enhanced  dissipation  entering  via  w,  rather 
than  thru  augmented  viscosity.  Unlike  the  Lumley[7]  scenario,  our  model 
does  not  require  the  incidence  of  a  coil-stretch  transition.  It  should  also  be 
mentioned  that  like  the  model  of  DeGennes  and  Tabor [8], [15],  our  model 
predicts  an  equipartition  of  turbulent  mechanical  and  elastic  energy.  It  does 
not,  however,  predict  a  ‘truncated’  cascade,  but  rather  a  modified,  viscoelas¬ 
tic  cascade.  Finally,  it  should  be  mentioned  that  our  model  fails  to  account 
for  the  observed  anisotropy  of  inertial  range  eddys  in  drag  reduced  flows  and 
for  the  relation  of  this  phenomenon  to  the  observed  decorrelation  of  stream- 
wise  and  cross-stream  velocity  fluctuations  and  thus  to  drag  reduction.  A 
possible  resolution  of  this  deficiency  may  come  from  the  proper  treatment 
of  the  small  scale  turbulence  dynamics  in  the  presence  of  a  strongly  sheared 
mean  flow. [25], [26] 
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4  CONCLUSIONS  AND  RECOMMENDA¬ 
TIONS 


4.1  Conclusions 

In  this  report,  we  have  surveyed  the  existing  experimental  and  theoretical 
understanding  of  drag  reduction  by  polymer  additives,  and  have  presented  a 
new  continuum  model  (polymer  hydrodynamics  —  PHD)  of  turbulent  poly¬ 
mer  solutions.  Several  natural  extensions  of  and  directions  for  further  devel¬ 
opment  of  this  model  are  apparent.  These  include: 

1.  the  application  of  standard  perturbative  closure  methods  (D.I.A.,  etc.)  [27] 
to  the  PHD  equations  to  develop  tractable  spectrum  evolution  equa¬ 
tions  for  PHD  turbulence. 

2.  investigation  of  PHD  turbulence  in  actual  bounded  shear  flow  configu¬ 
rations.  This  would  probably  develop  as  a  continuation  of  1.,  above. 

3.  the  generalization  of  the  PHD  model  to  incorporate  non-Hookean  elas¬ 
tic  eflects[16]. 

4.  an  investigation  of  the  coil-stretch  transition  in  the  context  of  a  dy¬ 
namically  relevant  chaotic  or  turbulent  flow. 

It  is  clear  that  the  study  of  drag  reduction  in  turbulent  polymer  hydro¬ 
dynamics  is  a  difScult  yet  fascinating  enterprise! 
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4.2  Program  Recommendations 


We  believe  that  the  time  is  right  for  an  experimental  and  theoretical 
program  to  understand  the  detailed  mechanics  of  polymer  drag  reduction. 
Such  a  program  would  consist  of  several  peurts: 

1.  Measurement  of  the  velocity  and  strain  rate  correlation  functions  and 
polymer  properties  (extension  and  orientation)  of  turbulent  polymer 
solution,  which  are  homogeneous  but  are  not  driven  at  an  outer  scale; 
homogeneous  turbulence  which  is  driven  in  a  Kolmogorov  cascade  from 
an  outer  scale;  and  pipe  wall  turbulence. 

2.  Measurement  of  drag,  turbulence,  and  velocity  structure  in  pipe  flows 
selectively  injected  with  polymers,  to  determine  which  effects  of  poly¬ 
mers  are  important:  creation  of  a  buffer  layer,  modiflcation  of  the  prop¬ 
erties  of  the  turbulence  in  the  fully  turbulent  fluid,  or  others. 

3.  An  extensive  program  of  numerical  simulation  of  turbulent  polymer  so¬ 
lutions  and  pipe  flows,  such  as  that  already  underway  at  the  La  Jolla 
Institute.  It  is  also  possible  to  make  diagnostic  measurements  (such  as 
of  the  Lagrangian  strain  rate  and  velocity  autocorrelation  functions) 
which  would  be  very  difflcult  in  the  laboratory.  The  great  flexibility 
and  versatility  of  these  methods,  complementing  real  laboratory  exper¬ 
iments,  makes  them  powerful  and  useful  tools. 

4.  Theory  to  understand  and  interpret  the  results  of  laboratory  and  nu¬ 
merical  experiments. 

On  account  of  the  difflculty,  complexity,  and  importance  of  polymer  drag 
reduction  phenomena  this  work  will  require  the  level  of  sophistication  in  ex¬ 
periment  and  theory  which  are  the  state-of-the-art  in  hydrodynamics,  plasma 
physics,  continuum  mechanics  and  computational  physics.  It  will  need  to 
draw  in  established  professionals  from  those  fields. 
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A  BASIC  PHENOMENOLOGY  OF  MHD 
TURBULENCE 

In  this  Appendix,  we  summarize  the  basic  phenomenology  of  MHD  tur¬ 
bulence  essential  to  our  discussion  of  polymer-hydrodynamics.  The  most  sig¬ 
nificant  feature  of  MHD  turbulence  is  the  modification  of  the  inertial  range 
spectrum  due  to  the  presence  of  ’large’  scale  magnetic  perturbations.  The 
effect  of  such  magnetic  perturbations  is  to  inhibit  inertial  range  transfer  by 
coupling  small  eddys  to  Alfven  waves  which  propagate  along  the  large  scale 
magnetic  perturbations.  This  transfer-inhibition  is  commonly  referred  to  as 
the  Alfven  effect  [20]. 


Here,  we  present  a  heuristic  derivation  of  the  Alfven  effect  by  exploring 
the  effect  of  a  uniform  magnetic  field  H  =  Boz  on  inertial  range  transfer. 
For  the  k  mode,  the  MHD  equations  are: 


~di 


^+(£•24  = 

-  ^(p  +  £.-&  +  ^/2)  (A-1) 


and 


^  +  (£  e4=b.|£*+(£  y£)^  (A-2) 


Here,  the  inertial  range  transfer  rate  for  the  eddy  of  scale  is  just  l/r*  = 
(K  •  2  which  we  estimate  using  standard  iterative  closure  methods,  i.e.: 


where: 


(A -3) 

*' 


Po 


=  i^'^kVk 
=  B..t.£g. 


(A-4) 


Here,  the  only  virtual  velocity-velocity  interactions  are  retained,  since  we 
only  seek  to  determine  the  effect  of  Bo  on  fluid  eddy-shearing.  Furthermore, 
for  simplicity,  a  single  parallel  wave  vector  is  assumed.  Here,  Awjt"  is  the 
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inverse  decorrelation  time  for  the  k  +  y  virtuail  mode.  Combining  Equations 
(A-2)  and  (A-3)  yields: 


1/  1 


+ 


(A -5) 


Here  =  Bo!  y/j^  is  the  Alfven  velocity.  Enforcing  locality  of  transfer  (i.e. 
At  fc')  then  yields: 


llTk^kVfJ 


(A -6) 


In  the  limit  — >  0  (i.e.  no  magnetization),  Equation  (A-6)  reduces  to  the 
familiar  relation  for  Kolmogorov  turbulence,  i.e.  l/r^  kVk-  In  the  strong 
magnetization  limit,  Equation  (A-6)  implies  that 

(A -7) 

It  is  apparent  that  a  strong  magnetic  field  (i.e.  Aw)  acts  to  reduce 

inertial  range  transfer.  Interpreting  Equation  (A-7)  iteratively  then  yields: 


KVa 


I/t*  ~  A;VJ 


ta 

^-(0) 


(A -8) 


Equation  (A-8)  states  that  a  strong  magnetic  field  inhibits  inertial  range 
transfer  by  increasing  eddy  lifetimes  by  a  factor  of 


The  result  of  Equation  (A-8)  is  appropriate  to  the  case  of  a  DC  magnetic 
field.  Since  large  scale  magnetic  fluctuations  will  behave  similarly,  and  since 
Alfvenic  interaction  is  non-local  in  k-space  (note  that  the  J.x.E.  force  cannot 
be  elinunated  by  a  Galilean  tra&sformationl),  we  can  generalize  Equation  (A- 
8)  by  arguing  the  correspondence  — ►  kVA,  where  Va  =  (5)rms/v^> 

the  Alfven  velocity  calculated  using  the  root-mean-square  fluctuating  field 
Bnosj  which  is  dominated  by  contributions  from  the  large  scale  magnetic 
perturbations.  In  that  case: 


1/t*  =  {kVuflkVA,  (A  -  9) 
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which  is  the  well-known  but  heretofore  unjustified  result  for  the  Alfvenically 
modified  eddy  turn-over  rate  appropriate  to  MHD  turbulence. 


The  notion  that  magnetic  fields  inhibit  inertial  range  interaction  is  fre¬ 
quently  expressed  by  the  alternative  statement  that  two  eddys  will  not  inter¬ 
act  unless  they  propagate  in  opposite  directions  along  neighboring  magnetic 
field  lines  (i.e.  see  Figure  3-5)[16].  This  conception  of  the  Alfven  effect  can  be 
manifested  in  the  simple  model  presented  by  relaxing  the  simplifying  assump¬ 
tion  of  a  single  A:,.  In  that  case,  A:,  — »  A:,-f  A:'  in  Equation  (A-6),  which  in  turn 
implies  that  energy  transfer  will  be  inhibited  unless  (A:,  -f  A:*  )VU  ^  Aw*#/.  Of 
course,  {kg +  k'^)Vji  — ♦  0  is  mathematically  equivalent  to  the  physical  picture 
that  two  Alfven  waves  propagate  in  opposite  directions.  It  is  also  interest¬ 
ing  to  note  that  the  ‘inhibition  factor’  ^1  +  {k“  k^Y j ^  suggests  an 
analogy  between  inertial  range  interaction  and  nonlinear  Landau  damping  in 
a  plasma[28],  with  fluid  eddys  as  the  analogue  of  particles  and  Au;j).»  as  the 
counterpart  to  the  spectrum-averaged  ballistic  frequency  (AA:)VJfc,  where  Vtu 
is  the  particle  thermal  velocity.  From  that  view-point,  the  inhibition  factor 
implies  that  wave-eddy  interaction  occurs  when  A:”  <  Aw*",  just  as  non¬ 

linear  Landau  damping  occurs  when  two  waves  beat  to  drive  a  ‘virtual’  mode 
which  resonates  with  the  particles,  i.e.  when  w  -|-  =  w"  <  AArVt/,.  Alter¬ 

natively,  just  as  spectral  transfer  due  to  nonlinear  Landau  damping  is  small 
in  the  ratio  {Ak)VthJuj,  inertial  range  transfer  in  MHD  turbulence  must  be 
small  in  the  ratio  AuicfkgVx’  In  the  corresponding  case  without  a  DC  field, 
AijjJkgVA  — »  V(^)/Va>  thus  recovering  Equation  (A-8). 


Having  derived  the  effective  eddy  turn-over  rate  for  MHD  turbulence,  it  is 
straightforward  to  determine  the  modified  inertial  range  spectrum.  Balancing 
the  dissipation  rate  e  (which  is  also  the  inertial  range  energy  flow-thru  rate) 
with  V{iyiT{()  (where  V{t)  and  T{t)  are  the  velocity  and  effective  turn-over 
rate  for  an  eddy  of  scale  /)  yields: 


r{i) 


vm^  M  iM 
e  T{i) 

iVA  ■ 


(A-IO) 
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Thus  V{i)  =  and  hence 


E{k)  =  (A -11) 


Thus,  the  Alfven  effect  changes  the  inertial  range  spectrum  from  E(k)  = 
^2/3^-5/3  E{k)  =  {tVAY^^k'^f^.  Note  that  since  rx(^)  <  t(£),  inertial 
range  mechanical  and  magnetic  energies  are  equal.  Similarly,  the  modified 
dissipation  scale  1^  can  be  determined  by  balancing  the  modified  eddy  turn¬ 
over  rate  with  the  viscous  dissipation  rate  where  u  is  the  kine¬ 

matic  viscosity.  This  implies: 

(A- 

so  that  using  the  relation  V{td)  —  {tVjYY^*^d*  yields 


!L^YML 

^  UVa 


U  =  (A  - 13) 


which  is  the  dissipation  scale  for  MHD  turbulence.  Note  that  the  assumption 
of  1/  (where  t)  is  the  resistivity)  is  implicit  in  this  analysis. 
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